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(N) 1 ABSTRACT 
+- > . 

O ■ We present the first complete ~55-671 pm spectral scan of a low-mass Class protostar (Serpens SMM1) taken with the PACS and SPIRE 

f™^ ■ spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12 CO (full ladder from 

W ■ /„=4-3 to 42-41 and EJk=491\ K), H 2 (up to 8, 8 -7o 7 and £ U /Jk=1036 K), OH (up to 2 n 1/2 7=7/2-5/2 and £ u /fc=618 K), 13 CO (up to 7„=16-15), 

Q\ \ HCNandHCO + (up to 7„ = 12-1 1). Bright [O i]63, 145 fim and weaker [C n] 158 and [C i]370, 609 ytim lines are also detected, but excited lines from 
chemically related species (NH3, CH + , CO + , OH + or H20 + ) are not. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here. 

' , '■ The - 10-37 yum spectrum has many fewer lines, but shows clear detections of [Nen], [Fen], [Sin] and [Si] fine structure lines, as well as weaker 

■ H 2 5(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (-54%), H 2 (-22%), [Oi] (-12%) and OH (-9%) 

T ^ ' emission. A multi-component radiative transfer model allowed us to approximately quantify the contribution of the three different temperature 

! components suggested by the 12 CO rotational ladder (r*°'w800 K, T™""*375 K and T™ ol *150 K). Gas densities m(H 2 )>5x10 6 cnT 3 are needed to 

i <~| reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope (warm and hot components) for which we derive upper 

Q_[ limit abundances of x(CO)<10 4 , x(H 2 O)<0.2xl0~ 5 and x(OH)<10~ 6 with respect to H 2 . The lower energy submm 12 CO and H 2 lines show 

I ' more extended emission that we associate with the cool entrained outflow gas. Fast dissociative /-shocks (v s > 60kms~') within an embedded 

, atomic jet, as well as lower velocity small-scale non-dissociative shocks (v s < 20kms~') are needed to explain both the atomic fine structure lines 

■4_> ■ and the hot CO and H 2 lines respectively. Observations also show the signature of UV radiation (weak [C 11] and [C 1] lines and high HCO + /HCN 

ffi " abundance ratios) and thus, most observed species likely arise in UV-irradiated shocks. Dissociative ./-shocks produced by a jet impacting the 
ambient material are the most probable origin of [Oi] and OH emission and of a significant fraction of the warm CO emission. In addition, H 2 
photodissociation in UV-irradiated non-dissociative shocks along the outflow cavity walls can also contribute to the [0 1] and OH emission. 
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1. Introduction us to detect a larger number of far-IR lines and to better con- 
C*~) ' strain their spatial origin. This is especially true for the detec- 
• ; In the earliest stages of evolution, low-mass protostars are deeply tion of high exc itation, optically thin lines that can help us to 
. embedded in dense envelopes of molecular gas and dust. During identify high temp erature components as well as to disentan- 
^ ■ collapse, conservation of angular momentum combined with in- gle and quantify the dominant heating mechanisms (mechanical 
fall along the magnetic field lines leads to the formation of a yfj radiative ). Early Herschel observations of a few low-mass 
dense rotating protoplanetary disk that drives the accretion pro- protostars for which partial or comp iete PACS spectra are avail- 
cess. At the same time, both mass and angular momentum are able show that the moleC ular emission of relatively excited 12 CO 
removed from the system by the onset of jets, collimated flows, (/<30)> j^q and QH lines is a common f eat ure (van Kempen et 
and the magnetic braking action (e.g., | Bontemps etalj [19961 aL 2010 for HH46 and Benedettini et al. 2012 for LI 157-B1). 
Ctf ; | Bachiller&Tafalla|1999|) . The resulting outflow produces a cav- This wam CO and h 2 emission was suggested to arise from 
ity in the natal envelope with walls subjected to energetic shocks the wallfj of an outflow . carved cavity in the envelope, which are 
and strong radiation fields from the protostar. These processes heated by UV pho tons a nd by non-d issociative C-type shocks 
heat the circumstellar gas and leave their signature in the pre- (yan Kempenetalj[2010; Vis ser etaljl2012h . The [Oi] and OH 
vailing chemistry (e.g., enhanced abundances of water vapour). emission observed in these sources, however, was proposed to 
Both the line and continuum emission from Young Stellar arise from dissociative /-shocks. Herczeg et al. (2012) detected 
Objects (YSOs) peak in the far-infrared (far-IR) and thus they higher-/ 12 CO lines (up to /=49-48) and highly excited H 2 
are robust diagnostics of the energeti c processes associated with lines in the NGC 1333 IRAS 4B outflow indicating the presence 
the first stages of star formation (e.g., |Giannini et aU|200l|), The of even ho tter gas. The very weak [O i]63 yum line emission in 
improved sensitivity and angular/spectral resolution of Herschel IRAS 4B outflow (the [0 1] 145 /mi line is not even detected) led 
spectrometers compared to previous far-IR observatories allows the se authors to conclude that the hot gas where H 2 dominates 

the gas cooling is heated by non-dissociative C-shocks shielded 



X 



* Herschel is an ESA space observatory with science instruments f r0 m UV radiation. Passive heating by the protostellar luminos- 
provided by European-led Principal Investigator consortia and with im- ity is also thought to contribute to the mid-/ 12 CO and 13 CO 
portant participation from NASA. 
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Fig. 1. PACS and SPIRE (lower panel) continuum-subtracted spectra of Serpens SMM1 (central spaxel) with H2O (blue), CO (red), 
OH (cyan), I3 CO (magenta), [0 1], [C 11] and [C 1] (green), and HCO + and HCN (black) lines labelled. The flux density scale is in Jy 
(assuming a point source). 



emission (IVisser et al.ll2012l lYildiz et al.ll2012l) . In NGC 1333 
IRAS 4A/4B, however, the 12 CO intensities and broad line- 
profiles of lower-/ transitions {J- 1 -0 up to 1 0-9) probe swept-u p 
or entrained shocked gas along the outflow (lYildiz et alJ l2012). 

In this paper we present the first complete far-IR and 
submm spectra of a Class protostar (Serpen s SMM1 or 
FIRS1), take n with the PACS (iPoglitsch et al.ll2010h and SPIRE 
(Griffin et al.ll2010[) spectrometer s on board the Herschel Space 
Observatory (Pilbrattet al. 2010). SMM1 is a low-mass proto- 



stellar system located in the Serpens co re dEiroa et al.l 2008) 



at a distance of c/=230+20pc (see also iDzib et al.l 20101 for 
an alternative measurement). It is the most massive low-mass 
Class source (M env ^16.1 M , r bo i^39K) of the WISH pro- 
gram sample {Water in Star-forming Regions with Herschel, 
Ivan Dishoeck et aD 1201 ll) . The Herschel spectra are comple- 
mented with mid-IR Spitzer Infrared Spectrograph (IRS) archive 
spectra of the emb edded protostar (photometry first presented by 
lEnoch et al.ll2009l) . 
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Fig. 2. Observed SED of Serpens SMM1 and modified black- 
body fit (red curve). Inset shows the full PACS array in flux den- 
sity units of 10~ 14 Wm^yUirr 1 . The abscissa is in linear scale 
from 55 to 190 /mi. The bright emission line seen in several po- 
sitions is [O i]63 /mi. Approximate outflow directions are shown 
with red and blue arrows. 



The Serpens cloud core was first mapped in the far-IR with 
the LWS spectrometer on board the Infrared Space Observatory 
(ISO). Even with a poor angular and spectral resolution of ~80" 
and R=A/ A/t-200, ISO-LWS observations of Serpens SMM1 re - 
vealed a rich spectrum of molecular lines ( Larsson et al.lj2 002). 
superposed onto a strong far-IR continuum Larsson et al.l2 000). 
In particular, 12 CO (up to /=21-20), H 2 (up to 4 4 i-3 30 ) and 
some OH lines were detected. Detailed modeling of the ISO 
emission concluded that those species trace the inner ~10 3 AU 
regions around the protostar, where gas temperatures are above 
300 K and densities above 10 6 cirT 3 . The comparison of the 
H2O, OH and CO lines with existing shock models suggested 
that shock heating along the SMM1 outflow is the most likely 
mechanism of molecular excitation. The details of these shocks, 
however, were less clear. 

Because of its far-IR luminosity and rich spectrum, 
Serpens SMM1 is an ideal Class source for studying the pro- 
tostellar environment and to characterize the shock- and UV- 
heating and chemistry with the much higher spatial and spec- 
tral capabilities provided by Herschel. In the following sections 
we present the data set, we identify the observed lines, present 
a 3 component non-LTE radiative transfer model and compare 
the observations with available shock models and observations 
of other YSOs. 

2. Observations and Data Reduction 

2.1. Herschel PACS and SPIRE Observations 

PACS spectra between ~55 and ~210/mi were obtained on 31 
October 2010 in the Range Spectroscopy mode. The PACS spec- 
trometer uses photoconductor detectors and provides 25 spec- 
tra over a 47"x47" field-of-view (FoV) resolved in 5x5 spa- 
tial pixels ("spaxels"), each with a size of ~9.4" on the sky. 
The resolving power varies between /?~1000 (at ~ 100 /mi in 
the Rl grating order) and ~5000 (at ~70/mi in the B3A order). 
The central spaxel was centered at the Serpens SMM1 proto- 
star (or 20 oo: 18 /, 29 m 49.8 i , <5 2 ooo: F15'20.5"). Observations were 
carried out in the "chop-nodded" mode with the largest chopper 
throw of 6 arcmin. The total observing time was ~3h (obser- 
vation IDs 1342207780 and 1342207781). The measured width 
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Fig. 3. Upper: SPIRE-FTS 12 CO 7=9-8 and 12-11 sparse maps. 
Lower. p-H 2 ln-Ooo and 2o 2 -ln maps respectively. The line 
surface brightness is on the 10~ 8 W m 2 sr 1 scale. The white 
contours represent the 50% line emission peak level. The 
FWHM beam is shown in each inset. Note the increase of line 
surface brightness in the more excited lines. 



of the spectrometer point spread function (PSF) is relatively con- 
stant for /l<100/mi (FWHM^spaxel size) but increases at longer 
wavelengths. In particular only ^40% (-70%) of a point source 
emission would fall in the central spaxel at ^190 /mi (A^6Qpm). 
Hence, extracting line surface brightness from a semi-extended 
source like Serpens SMM1 is not trivial. In this work we fol- 
lowed the method described by Herczeg et al. (2012) and Karska 
et al. (submitted). This involves multiplying all line fluxes mea- 
sured in the central spaxel with a /t-dependent correction curve 
calculated from the extended emission of the brightest CO and 
H 2 lines observed over the entire FoV. The correction factor is 
constant (~2) for /l<120//m and increases from ~2 to ~3-4 at 
~200 /mi. For the continuum and [0 1] and [C 11] lines, the indi- 
vidual line fluxes measured in all spaxels were added. Note that 
we checked that the reference spectra in the two off positions 
are free of [C 11] emission to confirm that the [C 11] lines towards 
Serpens SMM1 are real. Figure 1 shows the PACS spectrum in 
the central spaxel and Fig. 2 shows the entire PACS array from 
55 to 190 /mi (the bright emission line is [Oi]63 pm). 

SPIRE FTS spectra between -194 and -671 /mi (1545- 
447 GHz) were obtained on 27 March 2011. The SPIRE FTS 
uses two bolometer arrays covering the 194-3 13 /mi (Short 
Wavelength array, SSW) and 303-671 pm (Long Wavelength ar- 
ray, SLW) bands. The two arrays contain 19 (SLW) and 37 
(SSW) hexagonally packed detectors separated by ~2 beams 
(51" and 33" respectively). The unvignetted FoV is ~2'. The ob- 
servation was centered at Serpens SMM1 (IDs 1342216893) in 
the high spectral resolution mode (A/i=0.04 cm -1 ; fl~500-1000). 



3 



Goicoechea et al.: The complete spectrum of the Class protostar Serpens SMM1 





Ser SMM1: 63 fxm Ser SMM1: OH 65 




ARa. (arcsec) ARa. (arcsec) 

Fig. 4. PACS spectral maps of Serpens SMM1 in the o-H 2 3 3-2 12 (174.626/mi); [Cn] 157.741 yum; CO 7=17-16 (153.267 yum); 
both p-H 2 3 2 2-2 n (89.988/mi) and CO 7=29-28 (90.163 yum); [Oi]63.183//m and OH 2 n 3/2 7=9/2 + -7/2- (65.279 /an) lines (from 
top to bottom). The center of each spaxel corresponds to its offset position with respect to the protostar at a 2 ooo: 18 / '29'"49.8 S , 5 2 ooo: 
1°15'20.5". The Y-axis for each transition represents the normalized line flux (in the -0.2 to 1.2 range) with respect to the brightest 
line in the array. The X-axis represents the -550 to 550kms _1 velocity scale except for the OH 65.279 ^m line (-350 to 350kms _1 ) 
and the p-H 2 89.988//m and CO 90.163/vm lines (-1000 to lOOOkms 1 ). The vertical dashed lines show the position of the rest 
wavelengths. Blue and red arrows show the approximated outflow direction. The green box represents the central spaxel centered 
at the protostar position whereas the red and blue boxes represent postions in the red and blue outflow lobes respectively. Note that 
the PSF is relatively constant for /l<100jt/m (FWHM^spaxel size^9.4") but increases at longer wavelengths. 
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Fig. 5. (a) Spitzer IRS spectrum of the Class protostar Serpens 
SMM1 using the Hi Res modules and a 7"x7" aperture, (b) 
Continuum-subtracted spectra with main narrow lines identified 
in the LH module (A >19.5/mi) and (c) with the SH module 
(A<19.5pm). Note the lack of bright molecular line emission 
(H 2 Q, OH, etc.). 



The total integration time was ~9min. We corrected the 
spectrum to account for the extended continuum and molecu- 
lar emission of the source and the changing beam size with fre- 
quency across the SPIRE band. Taking account of the coupling 
of the source to the SPIRE-FTS beam (Wu et al. in prep.) and 
using the overlap region between the SSW and SLW bands as 
a reference (where the beam size differs by a factor of ~2) we 
calculate the FWHM size of the submm continuum source to 
be ~20". The SPIRE continuum fluxes were appropriately cor- 
rected and the unapodized spectrum was used to fit the line inten- 
sities with sine functions. Four representative mid-/ 12 CO and 
H2O integrated line intensity maps are shown in Figure 3. 

SPIRE and PACS data were processed using HIPE and then 
exported to GILDAS where basic line spectrum manipulations 
were carried out. Tables A.l and A. 2 summarize the detected 
lines and the line fluxes towards the protostar position (central 
spaxel; corrected for extended emission in the case of PACS ob- 
servations). The spatial distribution of several atomic and molec- 
ular lines observed in the PACS array is shown in Figure 4. 



2.2. Spitzer IRS Observations 

In order to complement the far-IR and submm Herschel spec- 
tra we have also analyzed mid-IR data towards the proto- 
star posit ion from the Infr ared Spectrograph (IRS) on board 
Spitzer dHouck et al.l 120041) in high spectral resolution mode 
(fl~600; 9.9-37.2//m). The Basic Calibrated Data (BCD) files 
for the Short High (SH) and Long High (LH) orders of IRS were 
retrieved from the Spitzer Heritage Archivefl (OBSID 34330, PI 

1 http://sha.ipac.caltech.edu/applications/Spitzer/SHA/ 



M. Enoch). Pipelin ing of the data was achieved using CUBISM 
(Smith et al.l20 07) to assemble the BCD files into spectral cubes. 
Bad and rogue pixels were removed using 3-sigma filtering. The 
spectra were extracted in a ~7"x7" aperture, thus close to the 
PACS spaxel size. The resulting Spitzer IRS spectra are shown 
in Figure 5. The photometry of these data was first presented by 
Enoch et al. (2009). Table A. 3 summarizes the detected lines and 
their fluxes. 



3. Results 

3. 1 . Far-Infrared and Submillimeter Lines 

Figure Q] shows the very rich PACS and SPIRE (lower panel) 
spectrum of the Serpens SMM1 protostar. More than 145 
lines have been detected, most of them rotationally excited 
lines from abundant molecules: 38 12 CO lines (up to 7=42- 
41 and £ u /fc=4971 K), 37 lines of both o-H 2 and p-H 2 
(up to 8i8-7i 7 and E U /A;=1036K), 16 OH lines (up to 2 U U2 
7=7/2-5/2 and £ u //fc=618K), 12 13 CO lines (up to 7=16-15 
and EJk=l 19 K) and several HCN and HCO + lines (up to 
7=12-11 and E U /&=283K). Weaker [Cn]158/mi and [Ci]370, 
609 pm lines are also detected. Excited lines from NH3, CH + , 
CO + , OH + or H 2 + are however not detected at the PACS 
and SPIRE spectral resolutions and sensitivity. The brightest 
line in the spectra is the [Oi]63/mi line (see also Figure 2) 
with a luminosity of L63-O.O14L . Assuming no extinction 
in the far-IR, we measure [Oi]63yum/[Cn]158jum-30 and 
[Oi]63jum/[Oi]145jum^lO line flux ratios integrated over the 
entire PACS array. 

Most of the far-IR continuum and line emission arises 
from the protostar position (central spaxel) with weaker, but 
detectable, contributions from adjacent outflow positions (see 
Figure 4). In particular, the high excitation lines of CO, H2O 
and OH detected below A<lQQpm seem to show a more com- 
pact distribution than the lower excitation lines (within the PACS 
PSF sampling caveats). The highest PACS spectral resolution 
(~90-60kms _1 ) is achieved in the ~60-70yum range. Although 
at this resolution all detected molecular lines are spectrally un- 
resolved (they have narrower intrinsic line-widths), the observed 
[Oi]63//m line is >30% broader than any close-by molecular 
line in the ~60-70yum range. This likely indicates that PACS 
marginally resolves the [Oi]63 pm line, being broader (or hav- 
ing higher velocity line-wing emission) than the excited far-IR 
CO, H2O and OH lines. In addition, the [Oi]63yum line-profile 
peak shifts in velocity from the protostar position (where the 
line is brightest) to the outflow red-lobe position where the line 
peak is redshifted by ~100kms~ 1 (more than a resolution ele- 
ment) and the molecular line emission is weaker (see red boxes 
in Figure 4). Note that [C 11] 158 /mi peaks at this red-lobe po- 
sition and the line peak is also redshifted by ~100kms _1 (al- 
though less than a resolution element in this wavelength range). 
Similar [Oi]63jt/m velocity-shifts are also seen in sources such 
as HH 46 and are associated w ith the emission from an atomic jet 
itself (Ivan Kempen et al] 1201 Oh . We refer to Karska et al. (sub- 
mitted) for a detailed comparison of the [O i]63 pm shifts in sev- 
eral outflows. 

Finally, the lower energy submm lines, from 12 CO and H2O 
in particular, show a more extended distribution than the far-IR 
lines, with some indication of outflow emission in the the S-E 
direction (see SPIRE maps in Figure 3). 

Table 1 summarizes the total line luminosities adding all ob- 
served lines between ~55 and 671 pm. Note that a few unidenti- 
fied lines (U) are present in the PACS spectrum. 
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Fig. 6. 12 CO and 13 CO rotational diagrams obtained from PACS 
and SPIRE spectra. The estimated error in T mt for 12 CO corre- 
sponds to different choices of J u cutoffs for the different compo- 
nents. The range of possible r rot for 13 CO reflects the standard 
error on the fitted values. 
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Fig. 7. H2O rotational diagram showing 0-H2O (filled circles) 
and P-H2O far-IR lines (open circles). Owing to subthermal ex- 
citation (T ml <s.Tk) and large opacities of most lines, the H2O dia- 
gram displays a larger scatter (standard error on the fitted values 
are shown) than CO. 



3.2. Mid-Infrared Lines and Extinction 

The mid-IR spectrum of Serpens SMM1 (Figure 5a) shows a 
much less rich spectrum than in the far-IR domain (probably 
due to severe dust extinction and lack of spectral resolution). 
Nevertheless, two weak H2 pure rotational lines (0-0 5(1) and 
5(2) transitions) and seven brighter atomic fine-structure lines 
of Ne + , Si + , S and Fe + (four transitions) are clearly detectecQ 
at the protostar positi on (Figures 5b and 5c). In contrast with 
NGC 1333-IRAS 4B dWatson et al. 1 120071) . no mid-IR OH and 
H2O lines are seen. The presence of [Ne n] and [Fe n] lines from 
energy levels above a few thousands Kelvin, together with the 
bright and velocity-shifted [0 1] line detected by Herschel, sug- 
gests the pr esence of fast dissociative sh o cks close to the pro- 
tostar {e.g., iHollenbach & McKed Il989b iNeufeld & Dalgarnol 
1989). Note that the same mid-IR lines ar e readily detected in 
Herbig-Haro objects and outflow lobes (e.g jNeufeld et all2 006; 
iMelnick et al.l 12008). As discussed in Sect. 4, the observed in- 
tensities of H2, H2O and CO rotational lines are more typical of 
slower velocity shocks. 

The large amount of gas and dust in embedded YSOs re- 
quires to apply extinction corrections to retrieve corrected mid- 
IR line luminosities. In this w ork we have taken gr ain optical 
properties from Laor & Draine dLaor & Draind [T993) and com- 
puted an extinction curve A(A)/Ay from IR to submm wave- 
lengths characterized by Rv-Av/E(B-V)— 5.5 (the ratio of visual 
extinction to reddening) and a conversion factor from extinc- 
tion to hydrogen column density of 1.4xl0 21 cirT 2 mag -1 . These 
values have been previously used to correct photome tric obser- 
vation s of embedded YSOs in the Serpens cloud (lEvans et al.l 
2009). Table 1 lists the uncorrected mid-IR line luminosities 
and those using an arbitrary correction of Ay=150. In the lat- 
ter case, even the [O i]63 fim line would be affected by extinc- 
tion (the de-reddened line flux would be ~2 times higher) and 
the [O i]63/145 /im line flux ratio would increase to ~20. 



2 First ionization potential of observed ionic species are: 21.56 (Ne), 
11.26 (C), 10.36 (S), 8.15 (Si) and 7.90eV (Fe). 



Such large intensity ratio s are predicted by disso- 
ciative shocks models ( e.g., [ Hollenbach & McKee 1989; 
Neufeld & Dalgarno 1989; Flower & Pineau Des Forets 2010) 
and thus provide a reasonable upper value to the extinction in 
the inner envelope. An independent estimate of the extinction is 
obtained from the ~9.7/im silicate grains absorption band seen 
in the Spitzer IRS low resolution spectra (shown in Enoc h et al.l 
120091) . In the diffuse ISM, the opti cal depth of the -9.7 fim broad 
feature is proportional to Ay (e.g., Whittet 2003). Although it is 
not clear th at the same correl ation holds in the dense protostel- 
lar medium (IChiar et al.l2007l) . we computed a lower limit to the 
extinction of A y>30 by fitting the ~9.7 jjm absorption band and 
assuming a single slab geometry. 

3.3. Dust Continuum Emission 

The dust continuum emission peaks in the far-IR domain at 
~ 100 /mi. This is consistent with the presence of a massive 
dusty envelope with relatively warm dust temperatures. Using 
a simple modified blackbody with a dust opacity varying as 
TU=Tioo(100/A)^, where T100 is the continuum opacity at 100 jum 
and B is the dust spectral index, we obtain a size of ~7" 
for the optically thick far-IR continuum source, with rioo— 2.5, 
B=\.l, Lfa,-_iR^26 L andadust temperature of 7^-33 K (Fig. 2). 
These parameters agree with more sophisticated radiative trans- 
fer models of the continuum emission detected with ISO/LWS, 
for which a far-IR source size of ~4" (~500AU in radius) 
and r,i=43K at the rioo=l equivalent radius were inferred 
(lLarsson et al.ll2000j). Thev are also cons istent with more re- 
cent SED models dKristense n et al .1 120 121) and with the com- 
pact (~5") emission revealed by mm interferometric observa- 
tions and thought to arise from the de nsest regions of the inner 
envelope and f rom outflow material (Hogerhe iide et al.1 1 1999: 
lEnoch et al.l2009t) . Note that owing to the large continuum opac- 
ity below 100/im, observations at these wavelengths probably 
do not trace the innermost regions, i.e., the developing cir- 
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Table 1. Observed and modelled luminosities. 
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<n) Fraction of total far-IR/submm line luminosity. 
*' From simple model: h=hot, w=warm, c=cool component. 
From Spitzer IRS observations (10-37 /mi). Upper limits include an 
arbitrary extinction correction of Av = 150mag. 



cumstellar disk with an approximate radius of ~100AU (e.g., 
iRodriguez etai1l2005h . 

4. Analysis 

4.1. 12 CO, 13 CO, H z O and OH Rotational Ladders 

Figure 6 shows all 12 CO and 13 CO detected lines in a rota- 
tional diagram. In this plot we assumed that the line emission 
arises from a source with a radius of 500 AU (see Sect. 3.3). 
Given the high densities in the inner envelope of protostars 
(>10 6 cm~ 3 ; see below), the rotational temperatures (r, ot ) de- 
rived from these plots are a good lower limit to 7\. The 12 CO 
diagram suggests the presence of 3 different T m components 
with 7: rot =622+30K (for J u <42), r rot =337±40K (/„<26) and 
r rot =103±15 K (7 U <14) respectively. The estimated error in T ml 
corresponds to different choices of J„ cutoffs for the different 
component^ In the following we shall refer to them as the hot, 
warm and cool components. Whether or not these T rot compo- 
nents are associated with 3 real physical components or with a 
more continuous temperature and mass distribution will be dis- 
cussed in Section 4. 1 . Here we note that the submm low-/ CO 
emiss ion shows broad line-profiles and a mo re extended distribu- 
tion (iDavis et alJl999HDionatos et alJ2010l) than the high-/ CO 
and H2O lines detected with PACS (more sharply peaked near 
the protostar). In addition, velocity-resolved observations with 
HIFI of high-/ CO lines up to /=16-15 show dif ferent profiles 
compared with low-/ CO lines dYildiz et al.ll2012l Kristensen et 
al., in prep). Specifically, for the 12 CO /=10-9 profile toward 
SMM1, Yildiz et al. (in prep.) find that about 1/3 of the inte- 
grated intensity is due to a narrow (FWHM of a few kms -1 ) 
component originating from the quiescent envelope and 2/3 to 
the broad outflow component. Thus, the cool gas component 
seen in the SPIRE submm maps is dominated by the entrained 
outflow gas. On the other hand, the CO /=16-15 profile of 
SMM1 is less broad and more similar to the excited H2O line 
profiles observed with HIFI (Kristensen et al., in prep.). Thus, 
the far-IR CO and H 2 lines detected by PACS clearly probe 
different physical components than the SPIRE data. 



3 Note that an extinction correction of Av/=150 will only imply 
slightly larger far-IR line fluxes (more than a factor 1.5 below ~80 yum). 
Slightly higher rotational temperatures would be inferred if Av/ = 150 
(~350 K and ~700 K for the warm and hot components respectively). 



Alternatively, Neufeld (2012) pointed out that the shape 
of the rotational diagrams do not, by themselves, necessar- 
ily require multiple temperature components. In particular, a 
single-TVot solution could be found for the 12 CO /„>14 lines in 
Serpens SMM1 if the gas were very hot (/ kin ^2500K) but had 
a low density (a few 10 4 ctrT 3 ). This scenario seems less likely, 
at least for the circumstellar gas in the vicinity of the protostar. 
Note that the gas density in the inner envelope of SMM1 is nec- 
essarily higher, as probed by our detection of high-/ lines from 
high dipole moment molecules such as HCN and HCO + (with 
very high critical densities). In addition, detailed SED models 
predict densities of «(H 2 )-4xl0 6 ctrT 3 at ~1000AU, whereas 
densities of the order of ~10 5 cm" 3 are only expected at much 
larger radii dKristensen et al.l 120121) . Therefore, it seems more 
likely that the different / rot slopes probe different physical com- 
ponents or the presence of a temperature gradient. 

A rotational diagram of the detected 13 CO lines provides 
Zrot=76+6K, thus lower than the T rot inferred for 12 CO in the 
cool component. The measured 12 CO/ 13 CO J-5-4 line intensity 
ratio is ^8, thus much lower than the typical 12 C/ 13 C isotopic 
ratio of ~60 lLanger & Penziasl (119901) . whereas the 12 CO/ 13 CO 
/=16-15 intensity ratio is ^55. These different ratios show that 
the submm 12 CO low- and mid-/ lines are optically thick, but 
the high-/ far-IR lines are optically thin. 

An equivalent plot of the observed far-IR H2O line intensi- 
ties in a rotational diagram gives / rot =l 36+27 K without much 
indication of multiple / lot components (Figure 7). Since H2O 
critical densities are much higher than those of CO, collisional 
thermalization only occurs at very high densities (>10 8 crrT 3 ). 
Therefore, the large scatter of the H2O rotational diagram is a 
consequence of subthermal excitation conditions (/rot^/k) and 
large H2O line opacities. We anticipate here that the opacity of 
most (but not all) far-IR H2O lines is very high (t»1). Compared 
to CO, the H2O rotational diagram is thus obviously less mean- 
ingful. A similar rotational diagram of the observed OH lines 
provides 7^=72+8 K (see also Wampfler et al., submitted), even 
lower than the inferred H2O rotational temperature. Like H2O, 
OH transitions also have high critical densities and line opacities 
(except the high-energy and the cross-ladder transitions). 

4.2. Simple Model of the Far-IR and Submm Line Emission 

Observations and models of the protostellar environment suggest 
that the observed emission can arise from different physical com- 
ponents, where different mechanisms dominate the gas heating 
and the prevailing chemistry: hot and warm gas from energetic 
shocks in the inner envelope (e.g., small scale shocks along the 
outflow cavity walls, bow shocks and working surfaces within 
an atomic jet); UV-illuminated gas in the cavity walls; cooler gas 
from the envelope passively heated by the protostar luminosit y 
and the entrained outflow gas (e.g., Ivan Dishoeck et all 1201 lb . 
Unfortunately, the lack of enough spectral and angular resolu- 
tion of our data does not allow us to provide a complete char- 
acterization of all the different possible components. However, 
the large number of detected lines, high excitation (£' U /A:>300 K) 
optically thin lines in particular, helps to determine the dominant 
contributions and the average physical conditions. 

With this purpose, we have carried out a non-LTE radia- 
tive transfer model of the central spaxel far-IR and submm 
spect rum (the protostar position) using a multi-molecule LVG 
code (Cernich arol2012l) . In this model we included three spheri- 
cal components suggested by the three temperature components 
seen in the 12 CO rotational diagram (i.e., hot, warm and cool 
components). In our model, CO is present in the three compo- 
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Table 2. Model components and source-averaged column densi- 
ties (assuming a line-width of 20kms _1 in all components). 



Comp. 


Radius 


r k (K) 


n(H 2 ) (cm- j ) 


Hot 


-500 AU 


800 


5xl0 6 


JV(cnT 2 ) 


12 CO(5xl0 lb ) 
13 CO(9xl0 14 ) 


H 2 O(2xl0 lb ) 




Warm 


-500 AU 


375 


5xl0 6 


Nicmr 1 ) 


12 CO(lxl0 18 ) 

OH(lxl0 16 ) 

HCN(2xl0 14 ) 

CO + (<2xl0 15 ) 

O(lxl0 18 ) 


1J CO(2xl0 lb ) 
HCO + (5xl0 14 ) 
CH + (<2xl0 14 ) 
NH 3 (<2xl0 15 ) 




Cool 


-4500 AU 


140 


2xlO i 


JV(crrr 2 ) 


12 CO(lxlO") 
13 CO(lxl0 16 ) 
HCN(lxl0 15 ) 
C + (lxl0 16 ) 


H 2 O(8xl0 15 ) 

HCO + (2xl0 14 ) 

C(lxl0 17 ) 





nents, but since a single rotational temperature roughly explains 
the observed far-IR H2O and OH lines, most of their emission in 
the PACS domain can be reproduced with only one component, 
the hot component for H2O and the warm component for OH 
(see next section for more details). Note that we do not model 
the extended quiescent envelope with the bulk of the (cold) 
mass seen in the mm continuum and narrow C ls O line emis- 
sion (Yildiz et al. in prep.). Foreground absorption and emission 
from the low density and low temperature extended envelope are 
thus not included in the model but they will basically influence 
narrow velocity ranges of the lowest energy-level lines. 

The latest available collisional rates were used {e.g., Daniel 
et al. 201 1 and references therein for H2O and Yang et al. 2010, 
extended by Neufeld 2012, for CO). Although only one line of o- 
H2 and one line of p-Ri are detected towards the Serpens SMM1 
protostar position, the large H2 0-0 S(2)/S{1) line intensity ratio 
suggests a H2 ortho-to-para (OTP) ratio lower than 3 in the gas 
probed by these low excitation H2 lines. Here we shall adopt an 
OTP ratio for H2 of 1 (the value we obtain from the two observed 
H2 lines assuming LTE and a rotational temperature of -800 K). 
Note that such low non-equilibrium H2 OTP values have been 
inferred, for example, in the h ot shocked gas towards HH54 and 
HH7-11 dNeufeld et al.ll2006h . The water vapour OTP ratio in 
the model, however, is let as a free parameter. 

The 0-H2O ground-state line observed with HIFI towards 
Serpens SMM1 shows a two component emission line profile 
with a medium component (Av^l5kms _1 ) thought to arise from 
small-scale shocks in the inner cavity walls and a broad compo- 
nent (Av ^40km s -1 ) from more extended shocked-gas along the 
outflow dKristensen et al.l [2012b . In addition, velocity -resolved 
observations of the OH 2 ni/2 7=3/2-1/2 line (-163. lpm) with 
HIFI show a broad line -profile with Av^20km s _1 (Wampfler, S. 
priv. comm.). For simplicity here we shall adopt a typical line- 
width of 20kms _1 in all modelled components. As in any LVG 
calculation, the model is more accurate for (effectively) optically 
thin lines (high-/ CO lines, excited and weak H2O and OH lines, 
etc.) than for very opaque lines {e.g., low excitation lines). 

4.2.1 . Model: Hot and Warm Component 

Following Sect. 3.3. and previous far-IR studies of this source, 
we have assumed a size of -4" (-500 AU in radius; e.g., 
lLarsson et aLlfeOOOl iRristense n et al. 2012) for the far-IR source 
(hot and warm components), a mixture of shocked and UV- 
illuminated gas (see Section 5). A good fit to the high-/ 12 CO 



emission is obtained for n(H 2 )~5xl0 6 cm 4 with 7\^800K and 
-375 K for the hot and warm components respectively. The 
approximate mass in these components is only M<0.03M o . 
Although this solution is not unique, this model satisfactorily re- 
produces not only the 12 CO high-/ lines but also the lines from 
other species (note that we searched for a combination of 7\ and 
n(H2) that reproduces all species simultaneously; see below). 

In order to constrain the density and place H2O and OH 
in a particular model component, we first run a grid of excita- 
tion models and investigated how particular line ratios change 
with n(H2) and 7\- In this comparison we choose faint, low- 
opacity lines arising from high energy levels that are observed 
with a similar PSF Figure A. 1 shows the variation of the 0-H2O 
3 3 o-3o3/6i5-3 5 (67.269/82.03 1/rni) and p-H 2 6 6-5i5/5 15 -4o4 
(83. 284/95. 627 yum) line ratios. The former ratio shows that the 
excited H2O lines arise from dense gas, whereas the latter re- 
quires high temperatures (>500K). The intersection of the two 
line ratios is consistent with densities higher than a few 10 6 cm -3 
and temperatures around -800 K. Note that for 7\<400K, 
the predicted emission from many excited H2O lines below 
-100 /im would be too faint. On the other hand, Figure A. 2 
shows a very clear dependence of the OH 65.279/96.363 pm line 
ratio with the gas temperature. The observed ratio does suggest 
that OH arises from gas at Tk<400K. Following the previous 
excitation analysis, and according to the higher observed rota- 
tional temperatures (r i ^°>/° t H ) and higher H2O energy levels 
detected, we placed H2O in the hot component and OH (together 
with [Oi]) in the warm component. Although this choice pro- 
vides a better fit to the data, we can not obviously conclude that 
the bulk of the water vapour and OH emission arise from differ- 
ent physical components. 

Far-IR radiative pumping is not included in our models and 
thus the derived H2O and OH column densities might be consid- 
ered as upper limits. Nevertheless, all H2O and OH lines are ob- 
served in emission and show r r ot>/d (where T& is the dust tem- 
perature inferred from the SED analysis, see Sect. 3.3), suggest- 
ing that owing to the high densities but moderate far-IR radiation 
field at >500AU, collisions dominate at first order. Note that 
high-mass protostars produce much stronger far-IR dust contin- 
uum fields (by orders of magnitude) and many far-IR OH and 
water vapour lines are often observed in absorption or show P- 
Cygni pr ofiles when far-IR pumping dominates the ex citation 
(see e.g., Goic oechea et a l. 2006; Cernicharo et al. 2006, for far- 
IR OH and H2O lines in Orion KL outflows). 

4.2.2. Model: Cool Gas Component 

In order to fit the more extended 12 CO low- and mid-/ emission 
and to match the !T ro t-100K component seen in the 12 CO rota- 
tional diagram (Figure 6), we included a third cool component 
- the entrained outflow gas. From the SPIRE 12 CO maps we in- 
fer an approximate radius of -20" (-4500 AU) for the emitting 
region. For this geometry, a density of n(H2)~2xl0 5 cirT 3 and 
a temperature of 7\-140K produ ces a good fit of the observed 
12 CO submm emission (see also lYildiz et al.ll2012l for velocity- 
resolved observations in NGC 1333 IRAS 4A/4B). This compo- 
nent is also needed to fit the lower excitation H2O submm lines. 

Low- and mid-/ 13 CO lines are optically thin and thus, 
in addition to the swept-up or entrained outflow gas, velocity- 
resolved submm 13 CO line profiles would carry information 
about other components such as the UV-heated gas and the qui- 
escent en velope, which may well dominate the 13 CO lower-/ 
emission (lYildiz et alJl2012h . Detailed models of the emission 
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Fig. 8. Synthetic non-LTE LVG spectrum of the Class protostar Serpens MM1 convolved to PACS (top three panels) and SPIRE 
(lower panel) spectral resolutions. Continuous curves show the line emission contribution from the "hot" (red), "warm" (green) and 
"cool" (blue) components discussed in the text. Vertical labels mark the wavelength position of CO (red), H2O (blue), OH (cyan), 
13 CO (magenta), HCO + and HCN (grey) rotational transitions. [Oi]63, 145 and [Ci]370, 609 yum fine structure lines are marked 
with green labels. 



from the extended and passively heated envelope (having most 
of the mass) predict 12 C Q rotational temperatures around ^30- 
60K (IVisseretal.ll2012L Harsono et al., in prep). 

To summarize, Table 2 shows the model parameters and 
Figure 8 shows the entire ~55 to 671 pm synthetic spectrum for 
the hot (red), warm (green) and cool (blue) components con- 
volved with the PACS and SPIRE resolutions. A comparison of 
the resultant synthetic spectrum (by simply adding the 3 compo- 
nent spectra) with the Herschel spectra is shown in Figure 9. 

4.2.3. Columns, Abundances and Validity of the Model 

Owing to the lack of angular resolution to determine the beam 
filling factors of the different physical components towards 
Serpens SMM1, the relative abundance ratios derived from the 
model and shown in Table 3 are a better diagnostic tool than 
the absolute column densities. Nevertheless, here we provide the 
source-averaged column densities (N) in the different compo- 
nents as well as the upper limits for several non-detected species 
(Table 2). In the hot component we find ~5xl0 16 cm~ 2 and 
~2xl0 16 cirT 2 column densities for 12 CO and H2O respectively. 
Note that a water vapour OTP ratio of ~3 provides the best fit 



to the observed water vapour lines and this is the value adopted 
in the models. Assuming that the gas in the warm component 
(7\-375K) covers a similar area, we obtain A^(OH)~10 16 cm -2 
and JV( 12 CO)~10 18 cirr 2 (a factor «20 larger than JV( 12 CO) in 
the hot component). 

The detected mid-IR H2 lines provide a lower limit to the 
H2 column density of A^h 2 ~10 22 cirT 2 . We use this column den- 
sity to provide an upper limit to the absolute abundances (with 
respect to H2) in the hot+warm components. In particular we 
obtain <10~ 4 , <0.2xl0~ 5 and <10- 6 for the CO, H 2 and OH 
abundances respectively. The inferred upper limit to the water 
vapour abundance is much higher than the ~10~ (8 ~ 9) value typi- 
cally found in cold interstellar clouds (e.g.. ICaselli et alj [2010) 
but is lower tha n the ^10~ 4 value often e xpected in the hot 
shocked gas (e.g jKaufman & Neu feld 1996). Note that the high 
source-averaged column density of atomic oxygen (~10 18 cm -2 ) 
suggests that a significant fraction of gas-phase oxygen reservoir 
in shocks is kept in atomic form. 

The column densities in the entrained outflow gas (cool com- 
ponent at 7\-140K) are given in Table 2. The absolute columns 
are more uncertain as they depend on the assumed spatial dis- 
tribution of the swept-up outflow gas. Note that in order to fit 
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Table 3. Selected abundance ratios in the modelled components. 



Species 


Abundance ratio 


Component 


H 2 0„/(CO„.+CO„) 


-0.02 


hot and warm 


OH,„/(CO, r +CO/,) 


-0.01 


hot and warm 


O lv /(CCv+CO„) 


-0.8 


hot and warm 


OH, t /H 2 0„ 


<0.5 


warm/hot 


co„/co,, 


-20 


warm/hot 


H 2 0„/CO„ 


-0.4 


hot 


ch ; ;/h 2 o„ 


<0.01 


hot 


0„/OH„. 


-100 


warm 


HCO+/HCN,,. 


-2.5 


warm 


HCO+/HCN, 


-0.2 


cool 


H 2 O c /CO c 


<0.08 


cool 



with observations suggests that this model captures the average 
physical conditions of the shocked gas near the protostar. The 
level of agreement is typically better than -30% (i.e., similar to 
the calibration uncertainty). The worst agreement occurs for sev- 
eral low-excitation optically thick lines that can arise from dif- 
ferent physical components because the non-local radiative cou- 
pling between different components is not treated in the LVG 
model. In addition, some lines su ch as the ~163^m OH lines 
(sensitive to far-IR pumping, se e lOffer & van Dishoeck 1992; 
iGoicoechea & Cernicharol [20021) are underestimated by >40%, 
suggesting that radiative pumping may play some role (see de- 
tailed OH models by Wampfler et al., submitted). 



the low-/ I3 CO lines in the cool component, we had to include a 
low 12 CO/ 13 CO column density ratio of ^10, confirming that the 
low- J 12 CO lines are optically thick and thus they do not probe 
the bulk of the material seen in the low-J CO isotopologue line 
emission (i.e., the massive and quiescent envelope). 

Our simple LVG model satisfactorily reproduces the abso- 
lute fluxes of most observed lines and does not predict lines that 
are not detected in the spectra (see Figure 9). Despite the fact 
that this model solution is obviously not unique, the agreement 



Regarding the 7\ and n(H2) conditions inferred in each com- 
ponent, Figure A. 3 in the Appendix shows that for a given gas 
density, we can distinguish temperature variations of -30% (es- 
pecially in high-/ CO and in some high excitation H2O and OH 
lines). Figure A. 4 shows that for a given gas temperature, density 
variations of a factor -2-3 can also be distinguished. Therefore, 
for the assumed geometry and physical conditions, the derived 
source-averaged column densities are accurate within a factor of 
-2. 
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5. Discussion 

5. 1 . Shocked-Gas Components 

Several fine-structure lines that probe the very hot atomic gas 
(a few thousand K) near the protostar are detected in the mid- 
IR. Besides, most tracers of the shock-heated molecular gas (a 
few hundred K) appear in the far-IR. At longer submm and mm 
wavelengths, extended emission from cool entrained outflow gas 
and from the cold massive envelope dominates. The relative in- 
tensities of the detected mid- and far-IR atomic and molecu- 
lar lines help to qualitatively constrain the nature of the main 
shocked-gas components in Serpens SMM1. 

The shock wave velocity, pre-shock density and the magnetic 
field strength determine most of the shocked gas properties. Fast 
shocks can destroy molecules and ionize atoms, whereas slower 
shocks heat the gas without destroying molecules. Depending on 
the evolution of the shock structure it is common to distinguish 
between /-type (or Jump) and C-type (or Continuous). More 
complicated, "mixed" non- stationary situations may also ex- 
ist (se e reviews by e.g., iDraine & McKedll993t [Walmsle v et al.l 
I2005I) . As we show below, our observations suggest the presence 
of both fast and slow shocks in Serpens SMM1. 



5.1.1. Fast Shocks 

The bright and velocity-shifted [Oi]63/mi emission, together 
with the detection of [Nen]12yum and very high energy [Fen] 
lines in Serpens SMM1 suggests the presence of fast dissocia- 
tive /-shocks related with the presence of an embedded atomic 
jet near the protostar. Note that [Ne n] and \Fe nl lines hav e been 
detected towards other YSOs (see e.g jLahuis et al.l2010l for the 
c2d Spitzer sample). Hollenbach & McKee (1989) presented de- 
tailed models for the fine-structure emission of atoms and ions 
in such dissociative /-type shocks. However, the chemistry of S, 
Fe, Si and related molecules (including gas-phase depletion, e.g., 
Nisi ni et al . 2005) were not included and thus the fine-structure 
absolute line intensity predictions are likely upper limits. In ad- 
dition, because the beam filling factors of the possible shock 
components are not known, line ratios provide a much better 
diagnostic than absolute intensities. According to these detailed 
models, the observed [Siii]35^m/[Fen]26/im <1 line intensity 
ratio (independently on the assumed extinction) provides a lower 
limit to the pre-shock density, «o=no(H)+2«o(H2), of >10 4 cm -3 . 
Besides, owing to the high ionization potential of Ne + , the 
[Nen]12/im/[Fen]26/im line intensity ratio increases sharply 
with the shock velocity (v s ). A lower limit of v s >60kms _1 is 
found from the observed [Ne n] 12 pm/[Fe n]26 pm >0. 1 line ra- 
tio (applying an extinction correction of Ay>30). 

Models of fast dissociative /-shocks predict that the gas 
is initially atomic and very hot (a few thousand K), but by 
the time that molecules reform, the gas cools to about 400 K. 
H2 formation provides the main heat source for this "tempera- 
ture p lateau" dHollenbach & McKedfl989T: iNeufeld & Dalgarnol 
1989). In such models, H2 is not so abundant and cooling by 
H2 and H2O can be less important than by [Oi], CO and OH. 
Therefore, in addition to the mid-IR fine-structure emission, a 
significant fraction of the [Oi], OH and the CO warm compo- 
nent emission can arise behind fast dissociative shocks triggered 
by a jet impacting the inner, dense evelope. 
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Fig. 10. Comparison of selected /7-H2O line intensity ratios ob- 
served with PACS and shock model predictions from Flower & 
Pineau Des Forets (2010). Red and blue curves represent /-type 
and C-type shock models respectively. The pre-shock gas den- 
sity is «o=2xl0 5 crrT 3 (continuous curves) and no=2xl0 4 cm~ 3 
(dotted curves). Low opacity lines from high energy rotational 
levels (£u/A:~300-650K) observed at similar wavelengths are se- 
lected. The upper panel shows the P-W2O 322-2n/6o6-5i5 line 
ratio (89.988/83.284jum) and the lower panel shows the 4 3 i- 
322/422-3o (56. 325/57. 637yum) line ratios. The grey horizontal 
areas show the observed ratios and their error margins. 

5.1.2. (UV-irradiated) Slow Shocks 

From the observed H2 0-0 S (2) weak line and the upper limit 
of the non-detected S (0) line, we derive a lower limit to the H2 
rotational temperature (T42) of ~700K (if a extinction correc- 
tion of Ay=30 is applied). This (and higher) temperatures of the 
mole cular gas are often inferred behind non-d issociative shocks 
(e. e.. iNeufeld et al.ll2006b iMelnick et al.ll2008h . 

Depending on the shock velocity, non-dissociative C-type 
shocks shielded from UV radiation can produce very high gas 
temperatures without destroying molecules (7\~400-3000K for 
v s ~10-40kms _1 in the models by Kaufman & Neufeld 1996). 
They also naturally produce high H2O/CO abundance ratios 
(close to 1) but predict low OH/H20<sl ratios (owing to neg- 
ligible H2O photodissociation but very efficient OH + H2 — > 
H2O + H reactions in the hot molecular gas). Therefore, a non- 
dissociative C-shock could be the origin of the observed H2 
lines, high-/ CO (/ u >30) and excited H 2 lines (£ U /A:>300K). 
Note that in our simple model we infer a high H2O/,/CO/j^0.4 
relative abundance when considering the hot component alone. 
On the other hand, if the hot and warm components inferred 
from the 12 CO rotational diagram only reflect the presence of 
a temperature gradient in the same physical component, then we 
would infer much lower H2O/,/(CO/,+CO, t )^0.02 abundances 
(see Table 3). Such low H2O abundances are difficult to recon- 
cile with non-dissociative C-shocks shielded from UV radiation 
and are more consistent with low-velocity /-shocks (see below). 

Planar-shock model predictions of absolute H2O line in- 
tensities depend again on the unknown beam filling factor of 
the shocked material. Figure 10 shows a comparison of se- 
lected H2O line intensity ratios with recent unidimensional mod- 
els of non-dissociative C-type and /-type shocks (the latter for 
v s <30kms _1 ) by Flower & Pineau Des Forets (2010). In order 
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to avoid using optically thick line diagnostics, different beam 
dilutions and extinction corrections, we selected low opacity p- 
H2O lines arising from high energy levels (E u / £~300-650 K) ob- 
served at similar wavelengths. The upper panel shows the 322- 
2n/6o6-5i5 line ratio (89. 988/83. 284/mi) and the lower panel 
shows the 43i-322/422-3i3 (56. 325/57. 637/mi) line ratio. The 
grey horizontal areas show the observed line intensity ratios 
and their error margins. From this comparison we conclude that 
low-velocity, non-dissociative shocks (v s <20km s -1 ) with a pre- 
shock density of «o^2xl0 5 cirT 3 are needed to reproduce the 
observed ratios. Figure 10 favors small-scale non-dissociative J- 
shocks, although it is difficult to determine any possible contri- 
bution of C-type shocks in the frame of current shock models 
and lack of spectral/angular resolution at /l<100//m. 

Besides, the very different shock velocities required by Ne + 
versus H2O lines suggests that both emissions arise from differ- 
ent locations. This will be the case if the Ne + emission arises 
from fast shocks within the jet itself, but the H2O (and H2) emis- 
sion arises in lower velocity shocks along the inner walls of the 
outflow cavity. Finally, note that the similar excitation conditions 
in the hot and dense gas (P/k- n Tk-8x10 9 KctrT 3 ), as well as 
the similar spatial distribution of the nearby CO J '=29-28 and 
/7-H2O 322-2n lines shown in Figure 4 suggests that the hot CO 
and the excited H2O lines arise in the same shock component. 

The UV radiation field probed by the [Cn]158//m line 
can also heat the gas along the outflow cavity walls (see e.g., 
IVisser et alj2012l) . UV photons will modify the chemistry of the 
shocked gas (e.g., by photodissociating H2O and enhancing the 
OH and O abundances; see Sect. 5.3), but in terms of gas heat- 
ing, they are likely not a dominant heating mechanism of the hot 
gas emission seen towards low-mass protostars (see Karska et al. 
submitted). All in all, the relative line intensities of the detected 
atomic and molecular species confirm the presence of both fast 
dissociative and lower velocity, non-dissociative shocks that are 
irradiated by UV radiation fields. 

5.2. Line Cooling 

Complete far-IR and submm spectral scans allow an unbiased 
determination of the line cooling in YSOs. Line photons follow- 
ing collisional excitation and escaping the region are responsi- 
ble of the gas cooling. Hence, the observed line luminosities 
in a broadband spectral scan provide a good measurement of 
the total gas cooling budget (see Table 1). In Serpens SMM1, 
~54% of the total line luminosity is due to 12 CO lines, fol- 
lowed by H z O (-22% and L(H 2 O)/L( 12 CO)^0.4), [Oi] (-12%) 
and OH (9%). The total far-IR and submm line luminosity 
iFiRSMM=i( 12 CO)+L(H 2 0)+L(Oi)+L(OH)+L( 13 CO)+L(C + )is 
£firsmm-O.12L , and the ratio L mo i/Lboi between the molecular 
line luminosity, L mo i=L( 12 CO)+L(H 2 0)+L(OH), and the bolo- 
metric luminosity is ^3.5x10 3 , all c onsistent with the ex pected 
emission of a Class source (see e.g jGiannini et al]20 01). Note 
that only if the extinction in the mid-IR is larger than Ay-200, 
then the intrinsic luminosity of the observed H2 lines will be 
higher than those of CO and H2O. 

In our model, CO dominates the line cooling (Table 1 ) with 
more than half of the total CO luminosity arising from the warm 
component (a mixture of shocks and UV-illuminated material) 
and with a warm/hot luminosity ratio of L(CO vv )/L(CO/,)-5. 
H2O is the second most important species, with a dominant 
contribution to the hot shocked gas cooling. Note that we pre- 
dict that >80% of the H2O line luminosity is radiated in the 
55//m</}<200yum range and only ~5% at shorter wavelengths 



(in agreement with the absence of strong H2O lines in the Spitzer 
observations). 

[Oi] and OH lines also contribute to the gas cooling. In 
particular, both the absolute OH luminosities and the observed 
[Oi]63//m/H2 0-0 5(1)»1 intens ity ratio are too bright for 
non-dissociative C-shock models dKaufman & Neufeldl ll996; 
Flower & Pineau Des Fbritsll2010h . As noted before, a signifi- 
cant fraction of the [0 1] and OH line emission likely contributes 
to the cooling behind a fast, dissociative /-type shock. 

5.3. Shock and UV-driven Chemistry 

The plethora of atomic fine -structure lines and to a lesser ex- 
tent, the relatively high OH/H2O<0.5 abundance ratio inferred 
in the warm+hot component, confirms the presence of strong 
dissociative shocks in the inner envelope. In addition, ow- 
ing to efficient H2O photodissociation, an enhanced UV radi- 
ation field illuminating the sh ocked gas can produce a high 
OH/H2O abundance ratios (cf. Goic oechea et al.1 1201 ll in the 
Orion Bar PDR). The UV radiation field near low-mass pro- 
tostars (roughly a diluted <1 0,000 K blackbody) is thought to 
be dominated by Lyar p hotons (10.2eV) a rising from accret- 
ing material onto th e star (|Bergin et all 2003) and from fast dis- 
sociative /-shocks (IHollenbach & McKedll979h such as those 
inferred in Serpens SMM1 (Section 5.1.1). Lya radiation can 
dissociate H2O (producing enhanced OH and O) but can not 
dissociate CO (~ll.leV) nor ionize sulfur (10.4eV) or carbon 
(1 1.3 eV) atoms. However, if significant C + and H2 are present, 
the ion-neutral drift in C-type shocks would significantly en- 
hance the CH + form ation rate compared to /-type shocks (see 
iFalgarone et all l2010i for the detection of CH + /=l-0 broad 
outflow wings in the shock associated with the DR21 mas- 
sive star forming region). The lack of CH + /=3-2 emission in 
Serpens SMM1 provides an upper limit for the CH + column 
density in the shocked gas (CH + /H2O<0.01) and will help to 
constrain UV- irradiated shocks models (Kaufman et al. in prep.; 
Lesaffre et al. in prep). The CH + /=3-2 line is detected in PDRs 
illuminated by mas sive OB stars whe re C + is the dominant car- 
bon species dGoicoechea et al.ll2.01 ll Nagy et al. in prep.) and, 
as in DR21, the UV radiation field contains photons with ener- 
gies higher than Ly a (up to 13.6eV). Excited lines from reactive 
ions such as CO + (that forms by reaction of C + with OH) and 
are related with the presence of strong UV radiation fields in the 
shocked gas, are not detected in the PACS spectrum despite pre- 
vious Jejitativ£_assignations of several high-energy far-IR lines 
(see lCeccarelli et alj[l997[ for IRAS 16293-2422) 

The spatial distribution of the [Cn]158//m line emission 
detected towards Serpens SMM1 (Figure 4) is similar to that 
of other species in the outflow (with brighter emission in the 
outflow position than towards the protostar itself). The detec- 
tion of faint C + emission indicates the presence of a relatively 
weak UV field (but able to ionize C atoms and dissociate CO) 
along the outflow. Besides, the increase of the HCO + /HCN 
abundance ratio in the warm temperature component (HCN is 
easily photodissociated while HCO + is abundant in the dense 
gas directly exposed to strong UV radiation) is another signa- 
ture of UV photons. Indeed, the mere detection of [Ci]370, 
609 pm lines at the poor spectral resolution of the SPIRE - 
FTS towards the protostar shows th at [Cil lines are signifi- 
cantly brighter than toward s HH46 (|van Kempen et al.l [2009) 
or NGC 1333 IRAS 4A/4B (lYildiz et al.ll2012l) . The presence of 
weak C + line emission and of Ly a radiation in protostars like 
Serpens SMM1 (although difficult to detect) suggests that UV- 
irradiated shocks are a common phenomenon in YSOs. 
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X-ray emission is a l so expected in low-mass YSOs 
(IStauber et al.ll2006l l2007t iFeigelsonl l20loh and they produce 
internally generated UV radiation fields after collisions of en- 
ergetic photoelectrons wit h H and H2 dMalonev et al.l Il996t 
iHollenbach & Gortil I2009I Although X-ray detections with 
Chandra have been reported t owards the nearby Se rpens SMM5, 
SMM6 and S68Nb protostars dWinston et al.ll2007l) . the strength 
of any X-ray emission from Serpens SMM1 is unknown, possi- 
bly because of the high column density to the embedded source. 

High sensitivity and velocity -resolved observations of CH + , 
CO + , SO + and other reactive ions related with the presence of 
ionized atoms (e.g., C + and S + ) will help us to characterize these 
UV-irradiated shocks in more detail. 

5.4. Comparison with Other Low-Mass Protostars 

Comparative spectroscopy of protostars in different stages of 
evolution allows us to identify the common and the more pe- 
culiar processes associated with the first stages of star forma- 
tion. Comparing the far-IR spectrum of the NGC 1333 IRAS 4B 
outflow to that of the Serpens SMM1 protostar, both show 
similar high H2O luminosities (L(H2O)~O.O3L ). However, 
Serpens SMM1 shows a factor ~50 stronger [Oi] luminosity 
and a lower L(H20)/L(CO) ratio. The weak [Oi] emission and 
high H 2 luminosity in IRAS 4B outflow (L(H 2 0)/L(CO)^l) 
has been inter preted as non-dissoci ative C-shocks shielded from 
UV radiation dHerczeg et al.ll20l2h . Indeed, C + is not detected 
in IRAS 4B outflow whereas it is detected in Serpens SMM1 
protostar and outflow (see Figure 4). On the other hand, the 
strong [Oi] and OH emission towards the Serpens SMM1 pro- 
tostar itself is more sim ilar to that of Class I source HH46 
(Ivan Kempen et alJl2010h where a L(OH)/L(H 2 O)^0.5 luminos- 
ity ratio has been inferred (vs. -0.4 in SMM1). One possibility 
is that a high-velocity jet impingin g on the dense inner env elope 
produces dissociative /-shocks (van Kempen et al. 2010) and 
enhanced [Oi] and OH emission compared to non-dissociative 
C-shocks. Although they could not distinguish the dominant sce- 
nario, either /-shocks or UV-irradiated C-shocks have been also 
proposed t o explain the O H emission seen in the high-mass YSO 
W3 IRS 5 dWampfleretal.ll20Tlh . 

The detection of [Nen], [Fen], [Sin] and [S 1] fine-structure 
lines towards Serpens SMM1 reinforces the scenario of both fast 
and slow shocks as well as UV radiation near the protostar, how- 
ever, it is difficult to extract the exact geometry of the different 
shock components in the circumstellar environment (e.g., Ne + 
versus H2O line emitting regions). Note that the presence of an 
embedded jet in the Class source L1448 was reported from 
the detection of [Fe 11] and [Si 11] lines by Dionatos et al. (2009), 
and the same lines have been dete cted towards t he nearby YSOs 
Serpens SMM3 and SMM4 dLahuis et al.ll2010t) . However, they 
did not d etect the fNen112.8/im line that requires fast sho ck 
velocities dLahuis et al.ll2007HBaldovin-Saavedra et al.ll2012l) or 
X-ray radiation (see Giidel et al. 2010, for Class II disk sources). 

Compared to HH46, CO lines with /„>30 and higher ex- 
citation H2O lines are detected in Serpens SMM1 (not neces- 
sarily due to different excitation conditions but maybe just be- 
cause lines are much brighter in SMM1). We have proposed 
that this hot CO and H2O emission arises from low velocity, 
non-dissociative shocks in the inner walls of the outflow cav- 
ity compressing the gas to very high thermal pressures (P/k- 
n 7k-10 9 ~ 10 KcrrT 3 ). Based on the large gas compression fac- 
tors and H2O line profiles seen toward several shock spots in 
bipolar outflows (far from the protostellar sources) Santangelo 
et al. (2012) and Tafalla et al. (in prep.) also conclude that cur- 



rent low velocity /-shocks models explain their observations 
better than stationary C-shocks models. New shock models with 
a more detailed geometrical description and including the effects 
of UV radiation are clearly needed to determine the exact nature 
of the shocks inferred in the dense circumstellar gas near proto- 
stars. 

The warm CO line emission (15</„<25) observed in 
Serpens SMM1 is a common feature in the far-IR spectrum of 
most protostars even if they are in differen t stages of evolution . 
This was previously observed by ISO (e.g., Gian nini et alJl200"Th 
and now by Herschel (e.g., Karska et al. submitted.; Green et al. 
in prep., Manoj et al., in prep.). This CO emission ubiquity nec- 
essarily means that a broad combination of shock velocities and 
densities can produce a similar warm CO spectrum in different 
low- and high-mass YSOs. 

The L(OH)/L(H 2 0), L(CO)/L([Oi]) and L(H 2 0)/L([Oi]) lu- 
minosity ratios in Serpens SMM1 have intermediate values 
between those in NGC 1333 IRAS 4B Class YSO and 
HH46 Class I YSO (with ratios closer to SMM1) suggesting 
that the different luminosity ratios in IRAS 4B outflow are either 
a consequence of an earlier stage of evolution, or simply because 
the IRAS 4B outflow is peculiar, with very high H2O abundances 
not affected by UV photodissociation. Indeed, one of the main 
conclusions of all these studies is that water vapour lines are 
a roboust probe of shocked gas in protostellar environments. In 
Serpens SMM1 we infer a total H 2 O/ 1 /(CO/ ] +CO w )^0.02 abun- 
dance ratio in the warm+hot components traced by our far-IR 
observations. Only if the hot component was an independent 
physical component, the H2O abundances would be much higher 
(H2O/,/CO/,-0.4). However, this value is higher than the typical 
H2O/CO abundance ratios inferred from HIFI observations of 
the ground-state 0-H2O 1 1 - 1 1 line (at ~557GHz) toward s a 
large sample of low-mass protostars dKristensen et al.lf2.Q12b . In 
these velocity -resolved observations, the H2O/CO abundance ra- 
tio increases from ~10~ 3 at low outflow velocities (<5 km s~') to 
-0.1 at high outflow velocities (>10 kms -1 ). Therefore, our in- 
ferred abundance ratio of H 2 O/ 1 /(CO; 1 +CO w )^0.02 from PACS 
observations seems a better description of the water vapour 
abundance in the protostellar shocked gas. Besides, it is also con- 
sistent with recent determinations from velocity-resolved HIFI 
obser vations of moderately excited H2O lines (EJk<215 K) (see 
e.g. JSantangelo et al.ll2012l forL1448 outflow). 

6. Summary and Conclusions 

We have presented the first complete far-IR and submm spec- 
trum of a Class protostar (Serpens SMM1) taken with Herschel 
in the -55-67 lyum window. The data are complemented with 
mid-IR spectra in the ~10-37ji/m window retrieved from the 
Spitzer archive and first discussed here. These spectroscopic 
observations span almost 2 decades in wavelength and allow us 
to unveil the most important heating and chemical processes 
associated with the first evolutionary stages of low-mass proto- 
stars. In particular, we obtained the following results: 

• Serpens SMM1 shows a very rich far-IR and submm emis- 
sion spectrum, with more than 145 lines, most of them rotation- 
ally excited lines of 12 CO (up to /=42-41), H 2 0, OH, 13 CO, 
HCN, HCO + and weaker emission from [C 11] 158 and [Ci]370, 
609 pm. [Oi]63jum is the brightest line. Approximately half 
of the total far-IR and submm line luminosity is provided by 
CO rotational emission, with important contribution from H2O 
(-22%), [Oi] (-12%) and OH (-9%). 
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• The mid-IR spectrum shows many fewer lines (probably 
due to severe dust extinction and lack of spectral resolution). 
Bright fine structure lines from Ne + , Fe + , Si + , and S, as well as 
weak H2 S(l) and S(2) pure rotational lines are detected. 

• The 12 CO rotational diagram suggests the pres- 
ence of 3 temperature components with -620+30 K, 
7^™~340±40K and r«*>'~ioo±15K. As predicted by SED 
models in the literature, the detection of H2O, OH, HCO + 
and HCN emission lines from very high critical density tran- 
sitions suggests that the density in the inner envelope is high 
(n(H 2 )>5xl0 6 cnT 3 ), and thus the CO rotational temperatures 
provide a good approximation to the gas temperature (r rot <7\). 

• A non-LTE, multi-component model allowed us to ap- 
proximately quantify the contribution of the different tempera- 
ture components (7^*800 K, T™""x375K and T™ ol *l50K) 
and to estimate relative chemical abundances. The detected mid- 
IR H2 lines provide a lower limit to H2 column density of 
Nh 2 > 1 22 cm~ 2 . We derive the following upper limit abundances 
with respect to H 2 in the hot+warm components: x(CO)<10~ 4 , 
x(H 2 O)<0.5xl0~ 5 and x(OH)<10~ 6 . The inferred water vapour 
abundance is higher than the value typically found in cold in- 
terstellar clouds but lower than that expected in hot shocked gas 
completely shielded from UV radiation 

• Excited CO, H2O and OH emission lines arising from high 
energy levels are detected (up to EJk-A91 1 K, EJk- 1 036 K and 
£u/A:=618K respectively). These species arise in the hot+warm 
gas (M<0.03 M Q ) that we associate with a mixture of shocks. 
The observed hot CO, H2O and H2 lines provide a lower limit to 
the gas temperature of ~700 K. The excited H2O line emission is 
consistent with detailed model predictions of low-velocity non- 
dissociative shocks (with v s <20kms _1 ). 

• Fast dissociative (and ionizing) shocks with velocities 
v s >60kms _1 and pre-shock densities >10 4 cirT 3 related with 
the presence of an embedded atomic jet are needed to explain 
the observed Ne + , Fe + , Si + and S fine-structure emission and 
also the bright and velocity-shifted [Oi]63//m line emission. 
Dissociative /-shocks produced by the jet impacting the ambient 
material are the most probable origin of the bright [0 1] and OH 
emission and of a significant fraction of the warm CO emission. 
In addition, water vapour photodissociation in UV-irradiated 
non-dissociative shocks along the outflow cavity walls can also 
contribute to the [0 1] and OH emission. 

Compared to other protostars, the large number of lines de- 
tected towards Serpens SMM1 reveals the great complexity of 
the protostellar environment. Both fast and slow shocks are 
needed to explain the presence of atomic fine structure lines and 
high excitation molecular lines. The spectra also show the sig- 
nature of UV radiation (C + and C towards the protostar and out- 
flow). Therefore, most species (including H 2 0, OH, O and trace 
molecules such as CH + ) arise in shocked gas illuminated by UV 
(even X-ray) radiation fields. Irradiated shocks are likely a very 
common phenomenon in YSOs. 
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Appendix A: Complementary Figures and Tables 



Goicoechea et al.: The complete spectrum of the Class protostar Serpens SMM1 



Table A.l. List of detected lines towards Serpens SMM1 and 
line fluxes from SPIRE and PACS central spaxels (corrected for 
extended emission in the case of PACS). The horizontal line 
near ~300 /im shows the transition between the SPIRE SLW and 
SSW detector arrays (beams between 42" and 17"). The line 
near ~200 fim marks the lines detected with PACS (spaxel size 
~9A"). Lines in the PACS ranges with detector leakage are not 
tabulated. Flux errors up to ~30%. Both the o-H 2 lot and p- 
H2O Ooo ground-state levels are at £] ow /fc=0K. 



Table A.2. Continuation of Table 1. 



Species 


Transition 


E«p 


A 


Line Flux 






(K) 




(W irT 2 ) 


CO 


19- 18 


1050.7 


137.196 


1.66E-15 


H 2 


5l4-5o5 


540.5 


134.935 


1.67E-16 


u 






132.982" 


1.13E-16 


u 






132.672° 


1.76E-16 


H 2 


4 23 -4i4 


397.9 


132.407 


3.95E-16 


CO 


20- 19 


1161.2 


130.369 


1.62E-15 


H 2 


4o4"3l3 


319.5 


125.353 


6.79E-16 


CO 


21 - 20 


1277.1 


124.193 


2.01E-15 


OH 


2 n 3/2 y=5/2 + -3/2- 


120.5 


119.441 


8.18E-16 


OH 


2 n 3/2 y=5/2--3/2 + 


120.8 


119.234 


8.75E-16 


CO 


22-21 


1398.6 


118.581 


1.29E-15 


H z O 


4l4"3o3 


289.3 


113.537 


7.39E-17 1 ' 


CO 


23 - 22 


1525.5 


113.458 


blended' 


CO 


24-23 


1657.9 


108.763 


1.08E-15 


H 2 


2 2 i-lio 


159.9 


108.073 


1.35E-16 


CO 


27-26 


2087.9 


96.773 


9.06E-16 


u 






96.633" 


6.14E-17 


OH 


2 n 1/2 - 2 n 3/2 y=3/2--5/2 + 


269.8 


96.363 


2.35E-16 


OH 


2 n, /2 - 2 n 3/2 y=3/2 + -5/2- 


270.2 


96.271 


2.35E-16 


H 2 


5l5"4o4 


470.0 


95.626 


5.16E-16 


H 2 


4 4 i-4 32 


668.1 


94.703 


1.76E-16 


CO 


28-27 


2242.2 


93.349 


7.55E-16 


u 






90.574" 


4.19E-17 


u 






90.434" 


4.64E-18 


CO 


29 - 28 


2401.9 


90.163 


6.57E-16 


H 2 


3 22 -2n 


296.8 


89.988 


8.29E-16 


CO 


30-29 


2567.0 


87.190 


7.03E-16 


u 






86.944" 


6.48E-17 


H 2 


7l6"7o7 


979.0 


84.766 


1.06E-16 


OH 


2 n 3/2 7=7/2--5/2 + 


290.5 


84.597 


9.29E-16 


CO 


31 - 30 


2737.6 


84.411 


1.25E-15 f 


OH 


2 n 3/2 y=7/2 + -5/2- 


291.2 


84.420 


blended 1 ' 


H 2 


6o6"5l5 


642.7 


83.283 


1.69E-16 


H 2 


6l6"5o5 


609.3 


82.031 


5.55E-16 


CO 


32-31 


2913.7 


81.806 


3.84E-16 


u 






81.013" 


7.03E-17 


CO 


33-32 


3095.1 


79.360 


2.49E-16 


OH 


2 n, /2 - 2 n 3/2 y=i/2 + -3/2- 


181.7 


79.179 


blended 1 


OH 


2 n, /2 - 2 n 3/2 y=i/2--3/2 + 


181.9 


79.116 


9.97E-16 1 


H 2 


4 23 -3i 2 


397.9 


78.741 


4.89E-17 


CO 


34-33 


3282.0 


77.059 


3.97E-16 


H 2 


3 2 i-2 12 


271.0 


75.380 


7.23E-16 


CO 


35-34 


3474.3 


74.890 


3.53E-16 


u 






69.840" 


5.75E-17 


CO 


38-37 


4083.5 


69.074 


1.57E-16 


CO 


39-38 


4297.3 


67.336 


1.28E-16 


H 2 


3o 3 -3o 3 


376.4 


67.269 


7.77E-17 


H 2 


3 3 i-2 2 o 


410.4 


67.089 


1.90E-16 


H 2 


3 3 o-2 2 i 


376.4 


66.437 


4.74E-16 


H 2 


7l6"6 2 5 


979.0 


66.092 


1.28E-16 


CO 


40-39 


4516.6 


65.686 


8.33E-17 


OH 


2 n 3/2 y=9/2 + -7/2- 


511.0 


65.279 


4.45E-16 


H 2 


6 2 5"5l4 


761.3 


65.165 


1.99E-16 


OH 


2 n 3/2 i=9/2--7/2 + 


512.1 


65.132 


4.11E-16 


CO 


41-40 


4741.2 


64.117 


1.23E-16 


H 2 


8o8"7l7 


1070.6 


63.457 


2.81E-17 


H 2 


8l8"7o7 


1036.5 


63.322 


4.00E-17 


O 




227.7 


63.183 


8.34E-15 


CO 


42-41 


4971.1 


62.624 


6.82E-17 


u 






62.003" 


9.05E-17 


H 2 


7 2 6"6l5 


1021.8 


59.986 


1.40E-16 


H 2 


4 32 -3 2 i 


516.1 


58.698 


2.59E-16 


H 2 


4 22 -3i 3 


454.4 


57.636 


2.64E-16 


H 2 


4 3 i-3 22 


552.3 


56.324 


1.05E-16 


U 






55.481" 


8.01E-17 



{a) Observed wavelength (otherwise rest wavelengths are tabulated) 



Species 


Transition 




A 


Line Flux 






(K) 


w 


(W irT 2 ) 


HCO+ 


5-4 


64.2 


672.326 


2.24E-16 


CO 


4-3 


55.3 


650.252 


4.11E-16 


C 


3 Pi- 3 Po 


23.6 


609.133 


8.69E-17 


HCO+ 


6-5 


89.9 


560.294 


3.16E-17 


13 CO 


5-4 


79.3 


544.161 


6.33E-17 


H 2 


lio-loi 


26.7 


538.289 


4.96E-17 


CO 


5 - 4 


83.0 


520.231 


5.13E-16 


HCO + 


7 - 6 


119.8 


480.275 


6.20E-17 


13 CO 


6-5 


111.1 


453.498 


1.11E-16 


CO 


6-5 


116.2 


433.556 


6.85E-16 


HCN 


8 - 7 


153.1 


422.912 


3.28E-17 


HCO+ 


8 - 7 


154.1 


420.263 


5.51E-17 


H 2 


2n-2o 2 


136.9 


398.643 


1.67E-16 


13 CO 


7 - 6 


148.1 


388.743 


1.42E-16 


HCO+ 


9-8 


192.6 


373.591 


4.10E-17 


CO 


7-6 


154.9 


371.650 


1.32E-15 


C 


3 P 2 - 3 P l 


62.5 


370.414 


1.42E-16 


13 CO 


8-7 


190.4 


340.181 


4.16E-17 


HCO+ 


10-9 


235.4 


336.255 


3.37E-17 


CO 


8-7 


199.1 


325.225 


1.68E-15 


HCN 


11 - 10 


280.7 


307.641 


5.55E-17 


I I f YV 


11 in 
1 1 - 1U 


OQO A 

zoz.4 


one n 1 
3Uj. / lz 


ftic n 
o.Vlr,-! / 


1 1 r\ 


1 

4)2" Ml 


lUU.o 


ana a^q 

DKJDADy 


A ACWt 1 
4.4UJ-.- 10 


H m 
CU 


n 
y - 


Id 1 .y 


1 r\o /lie 


1 enc 1 c 

l.oVb-lo 


tu 


o 
y - 


O/io ft 


oon ion 
zoV.lZU 


1 C./1C 1 c 




lz - 1 1 


DDD.O 


oon oco 
zoU.zj / 


O 1 OC 1 *7 


1 1 r\ 
I 1 2 v J 


"3 1. 

J 12-^03 


TIC T 

Zl J.Z 


07^. 1 QO. 
LID. LyD 


O.ZJJi- 10 




1 n n 
1U - y 


onn 


DID One 

z /z.zUj 


C AOTT 1 *7 

j.ozh,-! / 


1 1 ( \ 


1 n 

1 1 1 -Uno 


ddA 


OAO 07Q. 


/I /I |C 1 A 

4.41h,-lo 


tu 


in n 
1U - y 


3U4.2 


D£.r\ /in 
Z0U.Z4U 


no it 1 c 


1 1 r\ 

ti 2 U 


^21-^12 


Z/l.U 


0^7 7Qn 

Ld 1 . 1 y\) 


O /I QTh 1 A 
Z.4or> 10 


tU 


11 1 n 


LASl Q 

j'+o.y 


1A1 AQC\ 
Z4 / .47U 


1 ^ 1 P 1 A 


1 1 ( \ 


o o 
^20-^11 


1 CK Q 

Lyj.y 


DAT. Q70 
Z4j.y /Z 


O OOP 1 7 
Z.VZJi- 1 / 


t U 


11 1 n 
1 1 - 1U 




O^a ai a 

ZOO. 01 D 


O 1 71h 1 ^ 
Z. 1 /Ii-lZ? 


tu 


10 11 
1Z - 1 1 


A 1 O A 
41Z.4 


OOA GOQ 
ZZO.0V0 


7 A7P 1 7 
/ .0 /C.- 1 / 


t u 


19 11 
1Z - 1 1 


A~\ 1 3 
4 j 1. D 


1 (\ Q07 

zio.yz / 


Z.4 JJ_.- 1 J 


13 m 


1 "3 10 
YD - 1Z 


A Q. 1 n 
4ol .U 


ono /i7A 

ZUV.4 /O 


1 ^ 1 In 1 A 
1 .0 IJi-lO 


tu 


1 ^ 10 
YD ~ 1Z 


a 

jUo.O 


onn 070 

ZUU.Z /Z 


O AHTh 1 ^ 
Z.OUJi-lO 


13 co 


14- 13 


555.0 


194.546 


1.32E-16 


H 2 


4i 3 -4o4 


396.4 


187.110 


5.74E-17 


CO 


14 - 13 


581.0 


185.999 


1.55E-15 


13 co 


15 -14 


634.2 


181.608 


6.57E-17 


H 2 


2 2 ,-2i 2 


159.9 


180.488 


5.25E-16 


H 2 


2i 2 -loi 


80.1 


179.527 


9.72E-16 


H 2 


3o 3 -2i 2 


162.5 


174.626 


1.21E-15 


CO 


15 - 14 


663.9 


173.631 


2.20E-15 


13 CO 


16-15 


718.7 


170.290 


3.63E-17 


OH 


2 n 1/2 y=3/2--i/2 + 


269.8 


163.396 


1.92E-16 


OH 


2 n 1/2 i=3/2 + -l/2- 


270.2 


163.015 


2.76E-16 


CO 


16-15 


752.4 


162.812 


2.01E-15 


u 






161.808" 


3.49E-17 


c + 


2 p 2p 

"3/2- "1/2 


91.2 


157.741 


2.92E-16 


H 2 


322"3l3 


296.8 


156.194 


3.16E-16 


CO 


17 - 16 


846.3 


153.267 


2.15E-15 


O 


3 Po- 3 Pi 


326.6 


145.525 


8.38E-16 


CO 


18-17 


945.8 


144.784 


2.04E-15 


H 2 


3i 3 -2o 2 


204.7 


138.527 


9.49E-16 



(o) Observed wavelength (otherwise rest wavelengths are tabulated) 
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Table A.3. Mid-IR lines detected with Spitzer/IRS towards the 
Class protostar Serpens SMM1 and line fluxes (uncorrected for 
extinction) within a 7"x7" aperture. The horizontal line shows 
the transition between the LH and SH modules. Intensity errors 
up to -25%. 



OH 65.279/96.363 Intensity Ratio 



Species 


iransition 






A 


Line Flux 








(K) 


(41m) 


(WnT 2 ) 


Fe + 


b Ds/2- b Dj/2 




961 


35.349 


2.45E-17 


Si + 


2 Pil2- 2 P 1/2 




414 


34.815 


6.29E-17 


n 1 f„ 

p — JT12 


U U O \\J) J — z.- 


n 




28.219 


<5.92E-18 


Fe + 


6 Dj/2- 6 D 9 / 2 




554 


25.988 


8.12E-17 


S 


3 Pr 3 P2 




570 


25.249 


3.24E-17 


Fe + 


4 Z7 4 I? 
^5/2- -T7/2 




4087 


24.519 


4.46E-17 


Fe + 


4/7 4 17 
fj/2- rg/2 




3499 


17.936 


1.36E-17 


o-H 2 


0-0 5(1)7=3- 


-1 


1015 


17.035 


5.96E-19 


Ne + 


2 / J l/2- 2 / > 3/2 




1124 


12.814 


7.06E-18 


P-H 2 


0-0 5(2) 7=4- 


-2 


1682 


12.279 


2.42E-18 




o-H 3 3 30 -3 


33/6 




Intensity Ratio 



1400 



1200 
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log L0 n(H s ) (em ) 
p— H 3 S oe — 5 15 /5 L5 — Intensity Ratio 




0.4 



0.3 



0.2 



0.1 



log 10 n(H s ) (cm ) 



Fig.A.l. Grid of H2O LVG models for different gas tem- 
peratures and densities. Bottom: Contour levels of the p- 
H 2 6o6-5i 5 /5i 5 -4o4 (83. 284/95. 627/mi) line ratio for N(p- 
H2O)=10 16 cm -2 . Note the gas temperature dependence of this 
ratio. The black continuous curve shows the observed intensity 
line ratio of 0.33. Top: Contour levels of the 0-H2O 33o-3(b/6i5- 
3 05 (67.269/82.03 lpm) ratio for jV(o-H 2 O)=3xl0 16 cirT 2 . Note 
the gas density dependence of this ratio. The black continuous 
curve shows the observed intensity line ratio of 0.14 while the 
black dotted curve shows the observed /7-H2O 6o6-5is/5i5-4o4 
ratio. The intersection of both curves is marked with a star. 



2.5 







^^^^^^^ 






















Observed 






0.5 





- 2 



5 6 7 

log w n(H 2 ) (cnT 2 ) 

Fig.A.2. Grid of OH LVG models for different gas tempera- 
tures and densities. Contour levels of the OH 2 Hs/ 2 7=9/2 + - 
7/27 2 n 1/2 - 2 n 3 / 2 7=3/2--5/2 + (65.279/96.363/mi) line ratio for 
Af(OH)=10 16 cirT 2 . Note the gas temperature dependence of this 
ratio. The black continuous curves enclose the observed inten- 
sity line ratio with a ~ 15% of relative flux calibration error. 
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Fig. A.3. Best model discussed in the text (green curves) and a model where Tk is decreased by 30% in all components (blue curves). 
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Fig. A.4. Best model discussed in the text (green curves) and a model where n(H2) is decreased by a factor 3 in all components (blue 
curves). 
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